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We study the current statistics in normal diffusive conduc- 
tors in contact with a superconductor. Using an extension of 
the Keldysh Green's function method we are able to find the 
full distribution of charge transfers for all temperatures and 
voltages. For the non-Gaussian regime, we show that the 
equilibrium current fluctuations are enhanced by the pres- 
ence of the superconductor. We predict an enhancement of 
the nonequilibrium current noise for temperatures below and 
voltages of the order of the Thouless energy ETh ~ D/L^. 
Our calculation fully accounts for the proximity effect in the 
normal metal and agrees with experimental data. We demon- 
strate that the calculation of the full current statistics is in 
fact simpler than a concrete calculation of the noise. 

The electric current in conductors in general fluctuates. 
The full statistics of these fluctuations can be accessed in 
some cases. [Q However, the full statistics is not easily ac- 
cessible experimentally. The current experiments mainly 
concentrate on noise power measurements. This has be- 
come an important tool to extract information about the 
relevant charge transport mechanism. One can, for ex- 
ample, extract the effective charge of carriers responsible 
for the transport. Heterostructures with normal met- 
als (N) and superconductors (S) were intensively inves- 
tigated. Generally, the interest in these stems from the 
interplay between transport of doubly charged Cooper 
pairs and singly charged normal quasiparticles. The ra- 
tio between noise power in the superconducting state and 
in the normal state can roughly be interpreted as effec- 
tive charge. Experimental indications of a doubled shot 
noise due to Andreev reflection in diffusive wires with one 
superconducting lead have been reported in Refs. 
Other contributions report an enhancement of the cur- 
rent noise in SNS structures |^ with effective charges 
much larger that unity, possibly originating from multi- 
ple Andreev reflections. 

On the theoretical side first calculations of the noise 
[^,0 and the full statistics ||] for short contacts pre- 
dicted an enhanced shot noise with respect to the normal 
state value. For certain cases, such as tunnel junctions 
or disordered contacts, a doubled effective charge was 
predicted. The drawback of these calculations is the lim- 
itation to short contacts or small energies, since fully 
coherent propagation of electrons and holes is assumed 
and the energy dependence of the scattering amplitudes 



* Present address: Department of Physics, University of 
Basel, Klingelbergstr. 82, 4056 Basel, Switzerland 



has been disregarded. For the opposite regime, where co- 
herence of electrons and holes in the normal region plays 
no role, a modification of the Boltzmann-Langevin ap- 
proach has been put forward recently. ||^ Interestingly, it 
also gives a doubling of the shot noise in the incoherent 
regime. Similar limitations one finds in the available cal- 
culations of the noise of SNS structures. It is tempt- 
ing to "interpolate" between these two limits of coherent 
and incoherent propagation and to conclude, that noth- 
ing interesting happens in the intermediate regime. We 
will show below, that this is not the case. As analogy let 
us note the similarity to the linear conductance of a dif- 
fusive normal wire with one superconducting lead. Using 
the previously mentioned approaches one finds that the 
conductance has exactly the same value as in the normal 
state independent of temperature. Only a full calculation 
using the Keldysh Green's function technique revealed, 
that the conductance is significantly enhanced at energies 
of the order of the Thouless energy i?Th (l^ , displaying a 
reentrant behaviour. This behaviour of the conductance 
has been experimentally verified. At present no calcula- 
tion of current noise is available, which fully accounts for 
the peculiarities of the proximity effect. 

Another recent development is to study theoretically 
the full statistical properties of current fluctuations. The 
method of choice to do so is the so called full counting 
statistics, pioneered in Ref. Q. Introducing a counting 
fleld X that couples to the current operator, one can ac- 
cess the full current distribution. Derivatives of the cur- 
rent with respect to x immediately generate all moments 
of the distribution. Thus, one circumvents the cumber- 
some calculation of all the moments. 

In this Letter we introduce a theoretical method that 
allows to obtain the full current statistics for a wide class 
of SN-structures. We further present concrete results for 
the equilibrium current distribution and the nonequilib- 
rium current noise with full account of the proximity ef- 
fect. For simplicity, we restrict ourselves here to a dif- 
fusive wire mounted between a normal reservoir and a 
superconducting one. First we determine the equilib- 
rium current distribution. It turns out that the large 
current fluctuations are more probable in the supercon- 
ducting case. Second we calculate the temperature and 
voltage dependent current noise. We find that the noise 
is enhanced (similar to the conductance) at voltages of 
the order of and for temperatures below the Thouless 
energy. In the respective limits our results agree with 
previous calculation. [^|J^ 

Let us first present the theoretical framework. We use 
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a recently developed extension |^,|3[ of the Keldysh 
technique to compute the statistics of our proximity 
structure. There it was shown that the current statis- 
tics can be obtained by imposing the modified boundary 
condition on one reservoir's(L) Green's function 



Gl(x) 



(1) 



Here Gl is the standard Green's function of a reservoir 
and tk = o'lfs is a matrix in Keldysh(")-Nambu(~) space. 
Statistical properties are encoded in the dependence on 
the counting field x- The current probability distribution 
is found from 



P(/) 



dxe 



-S(x)-*X«/e 
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Here t denotes the time of observation. The action S{x) 
can be found from integration of the x-dependent current 



.e dSjx) 
't dx 



= Hx) = yJ ^^Tr [fKlix)] 



(3) 



The spectral matrix current in a diffusive wire is given 
by 



i{x) = -<^G{x)§^G{x) . 



(4) 



where a is the conductivity. This matrix current has to 
be found by defining an appropriate circuit. The manual 
on how to do this can be found in 1 14|. Note, that the only 
change in comparison to the calculation of the conduc- 
tance is the modified boundary condition (Q). All other 
relations defining the actual circuit remain unchanged, 
as long as they are respecting the full matrix structure. 

We stress that our approach is not in contradiction 
with the general scattering matrix approach of Ref . |^ . If 
one would know the electron-hole scattering amplitudes 
for the system under consideration, and would not disre- 
gard their energy dependence, one could obtain the same 
result. The characterization of these amplitudes would 
have to be performed along the lines of Ref. We 
also emphasize, that the calculation of the full current 
statistics can be done in a finer, simpler and more com- 
pact way than a separate calculation of it's second order 
perturbations series, i. e. the noise. 

Let us now specify our system. A diffusive metal is 
connected to a normal terminal on one end and to a 
superconducting terminal at the other end. Inside the 
mesoscopic wire the quasiclassical transport equations 
are obeyed. [|l6| In the normal metal they read 



D-^ {g{x,x)q^G{x,x) 



^ = [-iEf^,G{x,x)\ ■ (5) 



Here D is the diffusion constant and x the coordinate 
along the wire, which has a length L. It's conductance 



is Gn = crA/L (cross section A). At both ends bound- 
ary conditions to reservoirs have to be supplied. At the 
normal end with ideal connection the Green's function 
is continuous: (5(0, x) = Gl(x)- The other end is con- 
nected to a superconducting reservoir by a contact of 
negligible resistance, which leads to the boundary con- 
dition G(L,x) = Gr. a circuit representation of the 
system is depicted in the inset of Fig. |^. 

In a normal reservoir (which we will consider in the 
rest of the paper) Gl is given by 



Gt 



f:, K\ - fl-2f{E) 0. 
-fj = M l-2f{-E) ) ■ 



The distribution function at voltage V and temperature 
T is given by f{E) = iexp{{E + eV)/T) + A su- 

perconducting reservoir at zero voltage is described by 



R [R-A) tanh(£;/2r) \ 
Q A ) 



(7) 



Advanced and retarded Green's functions in (0) possess 
the structure R{A) = ffR.ATa + /n.An, /r,a = i^/{{E ± 
iSy — A^)^/^ and 5r,,a following from the normalization 
condition f^ + g"^ = 1 in a standard BCS-superconductor. 

Let us briefly comment on the numerical procedure 
of the solution. It is most convenient to solve the ma- 
trix equation (^) directly. For this purpose the diffusive 
wire is represented by a discrete set of n nodes, each 
represented by a Green's function Gk connected in se- 
ries by tunnel junctions of conductance g = {n + 1)Gn. 
I p^ The matrix current between two neighboring nodes 
is then given by Ik,k+i = f [Gk , Gk+i] ■ The right hand 
side of (|^) has a similar form and is represented as a 
'decoherence' terminal (with a Green's function Gdoc = 
—i2{E / nE'T\i)fz) connected to each node. The matrix 
current conservation for node k follows from discretizing 
Eq. (I) and reads [g{Gk+i + Gk-i) + Gdoc, Gk] = 0. The 
resulting set of equations for the nodes Green's functions 
is then solved by iteration. 

Because of its importance in the following discussion, 
we have depicted the differential resistance of the prox- 
imity wire in the inset of Fig.^ It displays the famous 
reentrance behaviour ||Tl[| , i.e. starting from the normal 
state resistance at zero voltage a minimum at energies of 
the order of several i?Tii occurs. Above the minimum it 
decays slowly to the normal state value as ^ (i?Th/£^)^^^. 
Therefore, the Thouless energy is the central quantity in 
the physics of the proximity effect. The disagreement of 
the theory and the experimental data from Kozhcvnikov 
et al. iQ may possibly result from heating effects, not 
accounted for in the theoretical calculation. 

As a first application, we study the distribution of cur- 
rent fluctuations in equilibrium. We therefore put V = Q 
and find the solution of the above equations for different 
values of x- We evaluate the integral over x in (|2|) in 
the saddle point approximation, we take x as complex 
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and expand the exponent in around maximum. The 
integral yields then P(/) « exp(5'(x) — x^t/e), which we 
plot implicitly as a function of /(x). To extract generic 
fluctuation properties of the proximity effect we set here 

A»r,£;Th. 

Results of this calculation are displayed in Fig. |l|. The 
current / is normalized by G-^T/e and lnP(/) is plotted 
in units of G-^Tt/e^. The solid line shows the distribu- 
tion in the normal state, which does not depend on the 
Thouless energy. In our units, this curve is consequently 
independent on temperature. In the superconducting 
state the Thouless energy does matter, and the distri- 
butions depend on the ratio Et\i/T. We observe that 
large fluctuations of the current in the superconducting 
case are enhanced in comparison to the normal case, and 
in both cases are enhanced in comparison to Gaussian 
noise. For comparison we plotted the Gaussian distribu- 
tion ^ cxp(— tJ^/4GNT) of the current measurements, 
owing to the fact that the conductance is the same in all 
cases. The differences between the normal and the su- 
perconducting state occur in the regime of non-Gaussian 
fluctuations. 




FIG. 1. Equilibrium current distribution. The case, 
when the superconducting terminal is normal (solid line) does 
not depend on the Thouless energy. In the superconducting 
state the distribution depends on Exh- At low temperatures 
(dashed curve) large fluctuations are enhanced in comparison 
to the normal state. This trend continues at temperatures 
above ETh (dotted curve). The deviations occur in the regime 
of non-Gaussian fluctuations. Note that the equilibrium noise 
is the same in all cases. The inset depicts our model system. 

Let us now turn to the noncquilibrium current noise 
as second application. Previous results can be summa- 
rized in two statements. At temperatures larger than 
an applied voltage the wire displays the usual thermal 
noise Si{T ^ eV) = AG^T in accordance with the 
fluctuation-dissipation theorem. Note that in a calcula- 
tion, which neglects the proximity effect, the conductance 



of the normal wire is equal to the normal state conduc- 
tance Gn independent of voltage or temperature. For 
large voltages eV ^ T, on the other hand, the shot noise 
Si{T <^ V) — (4/3)eGNV^ is doubled in comparison to 
the normal state value {2/3)eGj<iV. As discussed previ- 
ously these results do not depend on coherence between 
electrons and holes. Consequently the Thouless energy 
plays no role. However, this can not be true, since the 
conductance is enhanced by approximately 10% at ener- 
gies of the order of i?Tii- In the following we will show, 
that the proximity effect indeed changes the noncquilib- 
rium noise of the system. 

Using the method developed in 13|, we are able 



to take the coherence between electrons and holes fully 
into account. We have to solve the Usadel equation (|^) 
taking into account the boundary condition (^. Then, 
the current noise power is given by 



2ei 



dl{x) 



dx 



(8) 



There are two ways to attack this problem. One way to 
determine the noise is to expand the Green's functions 
and Eq. (|^) to first order in x ^nd thus to obtain an 
equation for the noise. The other way is to solve the 
full matrix equations, and perform the differentiation in 
(^) numerically. Below we will use the second way to 
find the noncquilibrium noise. Nevertheless, let us sketch 
the derivation of an equation for the noise. We define 
G(x,x)_ = Goix) - i{x/2)Gi{x) and /(a;,x) = Ux) - 
i{x/2)ii{x). As result we find 



h{x) = -a \G^[x) — G^{x) -f Gi(x) — Go(a;) 

-^hix)^{~iEf^,G^{x)\ . 
a ox 



(9) 



From this equation the generalization of the Boltzmann- 
Langevin equation to superconductors can in principle be 
derived. The boundary conditions at the reservoirs read 
Gi(0) = [-rK,GL] at the left end and Gi(L) = at the 
right end. Finally the noise is Si = —e J dETTf}^Ii{x). 
By taking the trace of Eq. ^ multiplied with tk it fol- 
lows that it does not matter, where the noise is evaluated, 
which is as it should be. 

We turn now to concrete results for the noise power. 
The influence of the proximity effect is most easily seen 
in the differential noise dS/dV. Numerical results for 
different temperatures are displayed in Fig. ||. The inset 
shows the differential conductance for the same parame- 
ters. The differential noise shows a remarkable enhance- 
ment at energies of the order of the Thouless energy. 
Following a linear increase at low voltage the differential 
noise overshoots the doubled normal differential noise, 
which is recovered at large voltages. The maximal dif- 
ferential noise occurs, if the voltage is of the order of 
« 4£'Th ■ The largest enhancement is found for T <C E'xh 
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and is roughly 10%. At higher temperature the differ- 
ential noise approaches the Boltzmann-Langevin result 
1^, shown as grey lines. At zero temperature the reen- 
trant behaviour makes the connection to the result ob- 
tained within random matrix theory. Q The nontrivial 
behaviour in the regime between these two approaches 
shows the importance of phase coherence. 




FIG. 2. Differential noise and conductance. Temperature 
and voltages are in units of -Eth. The main plot shows our 
results for the differential noise together with experimental 
data from Kozhevnikov et al. [Q| and the result of the Boltz- 
mann-Langevin approach (grey lines). The inset shows the 
differential resistance. Line styles and horiyontal axis corre- 
spond to the main plot, the small circles showing the corre- 
sponding experimental data of Kozhevnikov et al. . At 
zero temperature, both differential noise and resistance show 
a reentrant behaviour, featuring a distinct bump at voltages of 
the order of « 5i?Th. A finite temperature tends to smear this 
signature of the proximity effect, but an enhancement above 
the Boltzmann-Langevin result is still clearly visible. Com- 
parison with the experimental data show a good quali- 
tative agreement for T ~ 0.25i5Th, which corresponds to the 
experimental value. Note that our results have no adjustable 
parameter. 



We can compare our results with experimental data on 
the noise power obtained by Kozhevnikov et al.. [p|JT7| 
The experimental parameters correspond to T/E^h ~ 
0.2. Given that our approach contains no adjustable pa- 
rameters, the agreement is very good. A possible expla- 
nation for the difference between theory and experiment 
is that heating effects in the experiments may be impor- 
tant. These have been completely disregarded in the the- 
oretical calculation. This would also explain the smaller 
reentrance of the differential resistance seen in the exper- 
iment. Note, however, that the energy scale, at which the 
influence of the proximity effect is seen, is unchanged. 

In conclusion, we have developed a method to cal- 



culate statistical properties (with emphasis on the cur- 
rent noise power) of normal-metal-superconductor het- 
erostructures. The method is embedded in a matrix 
circuit theory jl^, which allows to find the full cur- 
rent statistics of a wide class of systems. We applied 
the method to a normal diffusive wire with one nor- 
mal and one superconducting reservoir. In equilibrium it 
turns out that the current fluctuation distribution differs 
from that of a purely normal system in the non-Gaussian 
regime. Large fluctuations of the current are enhanced 
by the proximity effect. For temperatures below i?Th 
we found that the nonequilibrium current noise shows 
a reentrant behaviour with a maximum for voltages of 
the order of -Exh- This is in qualitative agreement with 
experiment. 0| 
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